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We analyze a model system of fermions in a harmonic oscillator potential under the influence 
of a dissipative environment: The fermions are subject to a fluctuating force deriving from a bath 
of harmonic oscillators. This represents an extension of the well-known Caldeira-Leggett model to 
the case of many fermions. Using the method of bosonization, we calculate one- and two-particle 
Green's functions of the fermions. We discuss the relaxation of a single extra particle added above 
the Fermi sea, considering also dephasing of a particle added in a coherent superposition of states. 
The consequences of the separation of center-of-mass and relative motion, the Pauli principle, and 
the bath-induced effective interaction are discussed. Finally, we extend our analysis to a more 
generic coupling between system and bath, that results in complete thermalization of the system. 

PACS numbers: 03.65.Yz, 05.30.Fk, 71.10.Pm 



I. INTRODUCTION 

The study of quantum systems subject to dissipative 
environments is a topic which is both of fundamental im- 
portance in quantum mechanics and relevant for many 
applications requiring quantum-coherent dynamics. Fric- 
tion, energy relaxation, thermalization, destruction of 
quantum interference effects (decoherence) and the ir- 
reversibility of the measurement process are all examples 
of features that arise due to the system-environment cou- 
pling. The most important exactly solvable model in the 
theory of quantum-dissipative systems^ is the Caldeira- 
Leggett (CL) mode! 2,3 . It is the simplest quantum- 
mechanical model describing friction and fluctuations, 
and it has been used to analyze the quantum analogue of 
Brownian motion. The model consists of a single particle 
whose coordinate is coupled bilinearly to the coordinates 
of a bath of harmonic oscillators. Thus, it can be solved 
exactly as long as the particle either moves freely (apart 
from the system-bath coupling) or inside a harmonic os- 
cillator potential. 

However, in several important applications of these 
concepts we are not dealing with a single particle but 
rather with a many-particle problem from the outset. 
This applies in particular to solid state physics, where 
the dephasing of electronic motion often determines the 
temperature-dependence of quantum-mechanical inter- 
ference effects such as Aharonov-Bohm interference, uni- 
versal conductance fluctuations, or weak localization. 
Usually, the fermionic nature of the electrons makes it 
difficult to apply the insights gained from single-particle 
calculations, since the Pauli principle may play an im- 
portant role in relaxation processes. 

In this paper, we present our results for a natural ex- 
tension of the CL model to a case with many fermions. 
We consider a number of non-interacting fermions mov- 
ing inside a harmonic oscillator potential, subject to a 



Figure 1: Left: Coupling the fermions in the oscillator to 
the fluctuating quantum force F leads both to relaxation of 
excited fermions (decay rates depending on the bath spectrum 

(^FFj ), as well as some smearing of the level occupation 

even at T = 0. Right: The bosonized version is equivalent 
to chiral fermions moving on a ring, subject to a transverse 
fluctuating force. 



fluctuating quantum force deriving from an oscillator 
bath (Fig. QJ. In the single-particle case, this is just the 
CL model, describing a quantum damped harmonic os- 
cillator. Thus, we are studying a "Fermi sea in a damped 
harmonic oscillator" . Our main motivation is to analyze 
an exactly solvable model problem featuring dephasing 
and relaxation in presence of the Pauli principle. Never- 
theless, the model may also become helpful for the study 
of cold fermionic atom clouds in quasi one-dimensional 
harmonic atom traps, where fluctuations of the external 
trapping potential might be described by a homogeneous 
fluctuating force acting on the cloud. Alternatively, it 
could be an approximate description of a parabolic (quasi 
one-dimensional) quantum dot subject to the Nyquist 
noise of a fluctuating electric field from nearby gates. 

There are only a few previous theoretical works study- 
ing dephasing and relaxation of fermions in the context of 
quantum-dissipative system (excluding Boltzmann-type 
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kinetic equations or simple perturbation theory). To the 
best of our knowledge, a many-fermion generalization of 
the CL model was first suggested in Ref. 0, where free 
fermions coupled to independent oscillator baths were 
studied. As this breaks the indistinguishability of the 
fermions, it is not clear to which extent this model may be 
realized physically, and the influence of Pauli-blocking on 
relaxation processes could not be studied. A discussion 
of Luttinger liquids coupled to dissipative environments 
has been provided in Ref. where the emphasis was 
on general features rather than the actual evaluation of 
Green's functions. Some aspects of our model are similar 
to the features found in a study& of interacting fermions 
in a parabolic trap without coupling to a bath. Another 
extension of the CL model to many fermions^ involved 
particles on a ring subject to a quantum force indepen- 
dent of the position on the circumference. However, this 
did not reveal any influence of the Pauli principle, as such 
a coupling does not lead to transitions between different 
momentum states (Pauli blocking becomes relevant only 
for tunneling from external Fermi systems). Coupling 
a heat bath to indistinguishable fermions has also been 
employed as a tool for quantum molecular dynamics 8 . 
For a recent detailed perturbative study of Nyquist noise 
leading to decoherence of electrons in quantum dots and 
wires, for realistic gate geometries, we refer the reader 
to Refs. U The Feynman- Vernon influence functional, 
widely used in the analysis of a single dissipative parti- 
cle, cannot be applied directly to the many-fermion sit- 
uation. In the theory of weak-localization, there have 
been phenomenological prescriptions to incorporate the 
effects of the Pauli principle in an influence functional 
approach^. Recently, a formally exact generalization to 
the many-fermion cas o 11 ' 12 has been derived, although 
the evaluation of the resulting path-integral remains dif- 
ficult. 



II. THE MODEL 

We consider a system of N identical fermions (non- 
interacting and spinless) confined in a one-dimensional 
harmonic oscillator potential. They are subject to a fluc- 
tuating force F that by necessity acts on each fermion 
with equal strength, i.e. it couples to the center-of-mass 
coordinate of the system of fermions (see Fig. [JJ, like 
oc F^2j Xj. This force derives from a bath of oscillators. 
In second quantiziation, the Hamiltonian reads: 



H = u! ^ ncic n + H B + 

oo 

^ Vn + l(cjj +1 c„ + h.c.) 



F 



(1) 



The c n are fermion annihilation operators, and the oscil- 
lation frequency of a fermion of mass m in the parabolic 
potential is luo (we have set h — 1). Note that uq already 
contains a counterterm that depends on the coupling to 
the bath, see Eq. (JSJ) below. The bath Hamiltonian is 
given by: 



N B 

£ 



(2) 



Here the bath oscillator masses have been chosen to be 
equal to that of the fermions, without any loss of gen- 
erality, in order to streamline a few expressions derived 
in the next section. The force F is given as a sum over 
the bath normal coordinates Qj (with a prefactor g of 
dimensions energy over length squared): 



The work is organized as follows: After introducing the 
Hamiltonian in Sec. [H] we apply a certain approxima- 
tion in order to rewrite and solve the Hamiltonian using 
the method of bosonization, in Sec. IIIII This forms the 
basis for our evaluation of the fermions' single-particle 
Green's functions, described in Sees, and leading 
to the exact general expression in Eqs. (|38|l . I|39|l . In 
Section IvTI we analyze in more detail our approximation 
of 'large' particle numbers needed for bosonization. In 
Section TVIII we evaluate the two-particle Green's func- 
tion (exact result in Eqs. (|B7(l -i|B9 |) L which enables us 
to discuss the decay of populations after adding an extra 
single particle above the Fermi well as dephasing 

of a coherent superposition of states. Finally, we extend 
our model to a more generic type of coupling between 
fermions and bath (Sec. IVIII(I . Technical details are rel- 
egated to the appendices. 



N B 



F 



(3) 



Its spectrum 



+ 00 



F{t)F)e +lult dt 



(4) 



is still arbitrary and depends on the distribution of bath 
oscillator frequencies. At T = 0, it is given by: 



(ff) =— V SUa-Slt). (5) 



The special case of an Ohmic bath, as it is used in the 
theory of Quantum Brownian motion 2 , is defined by 



Some of these results have already been presented in a 
brief version^. 



T=0 



(PfY U = ^-ujd{uj c - oj)9{lu) , (6) 
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where rj = 7717 is the coefficient entering the friction force 
— T]v acting on a single particle of mass m, with 7 the 
corresponding damping rate. 

The Hamiltonian (JJJ can be derived from the follow- 
ing form, where the fcrmions are treated without second 
quantization and the translational invariance of the cou- 
pling between fermions and bath particles is apparent: 



N 



2m 



N N E 



Jul 

2 

N 



Xf + 



Pi , mW 00~2 



E 



1=1 3=1 



2m 



(7) 



Here the couplings gj are given in terms of the parameter 
g as <7j = g 2 N 2 /(2mf2|ATs), and the rescaled bath co- 
ordinates are Qj = Qj (mil 2 y/Ns /(gN)). Note that the 
frequency loq introduced above already contains the coun- 
terterm that arises from the xf -terms in (JJJ and which is 
essential to prevent the effective oscillator potential from 
becoming unstable for larger coupling strengths g. In 
terms of the bare oscillator frequency loqq, it is given by: 



^00 



2— / dijj 
m Jo 



I - - \ T=0 



(8) 



III. 



SOLUTION BY BOSONIZATION AND 
DIAGONALISATION 



Since the fluctuating force acts only on the center-of- 
mass (cm.) coordinate of the particles and, in the case 
of the harmonic oscillator potential, the cm. motion is 
independent of the relative motion, the model defined 
above is, in principle, exactly solvable in a straightfor- 
ward manner. The solution can be carried out by finding 
the classical eigenfrequencies and eigenvectors of the to- 
tal system of N + Nb coupled oscillators, setting up the 
quantum-mechanical wave functions in the total Hilbert 
space and performing the antisymmetrization with re- 
spect to the fermion coordinates. Afterwards, any desired 
observable (reduced density matrix, occupation numbers, 
fermion Green's functions etc.) may be calculated in 
principle. However, due to the antisymmetrization and 
the appearance of the oscillator eigenfunctions in the in- 
termediate steps of the calculation, this procedure gets 
extremely cumbersome, so we prefer a different route 
which is approximately valid for large fermion numbers 
N. In Section IvTl we will analyze some aspects of the 
small N case and compare with the bosonization results. 

We assume the number N of fermions in the oscillator 
potential to be so large that the lowest-lying oscillator 
states are always occupied, for any given many-particle 
state that becomes relevant in the calculation (at the 
given interaction strength and temperatures). In other 



words, the excitations in the fermion system, induced 
by the bath (and temperature), are confined to the re- 
gion near the Fermi surface. Then we may employ the 
method of bosonization where the energy of the fermions 
is rewritten as a sum over boson modes (i.e. sound waves 
in the fermion system). This is possible since the ener- 
gies of the oscillator levels increase linearly with quantum 
number, just as the kinetic energy of electrons in the Lut- 
tinger model of interacting electrons in one dimension. 
The same procedure has been applied to describe inter- 
acting fcrmionic atoms in a parabolic one-dimensional 
trap, for certain exactly solvable model interactions*. 
For recent pedagogical reviews on Luttinger liquids and 
bosonization, see Ref. 0. We note that this equal spac- 
ing of energy levels is approximately present for any po- 
tential at high excitation energies, where a semiclassical 
description of the single-particle dynamics becomes valid. 
However, the applicability of the following description to 
such a situation also depends on the structure of matrix 
elements of the fluctuating potential in the single-particle 
eigenbasis. For the harmonic oscillator considered here, 
only adjacent levels are connected by the position opera- 
tor x (describing the potential of the homogeneous force) . 

We introduce (approximate) boson operators, which 
destroy particle-hole excitations (q > 1): 



^ 00 

~~f=j E C n C n+q 
V 9 „=0 



(9) 



One may check that they fulfill the usual boson commu- 
tation relations, up to terms involving levels near n = 
that vanish when acting on the many-particle states oc- 
curing under our assumptions. Alternatively, one may 
make these relations hold exactly by redefining the orig- 
inal model to incorporate an infinite number of artificial 
single-particle levels of negative energy, as is done in the 
Luttinger model. 

Then the central result of bosonization may be applied, 
i.e. the fermion energy (bilinear ) may be written as 
a bilinear expression in b q . This is possible only due to 
the linear dependence of energies on the quantum number 



loo nclc r = 

n=Q 



" = Ua E q H^1 + E N 



(10) 



9=1 



Here = ujqN(N — l)/2 is the total energy of the 
7V-fermion noninteracting ground state. We keep N as 
an operator at this point since we will be interested in 
calculating Green's functions where the particle number 
changes. 

Under the same assumption of large N 3> 1 we get: 



^ Vn + l(c]* l+1 c„ + h.c. 



71=0 
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Again, this has to be understood as an approximate op- 
erator identity which is valid when applied to the many- 
particle states we are interested in, where n « N. We 
approximate N to be a number instead of an operator 
in this formula. This expression shows that the fluctu- 
ating force only couples to the lowest boson mode (with 
<7=1), corresponding to the cm. motion. Therefore, the 
Hamiltonian now has become (within the approximations 
described above): 



with 



H w ujq ^2 qb\b q 



N 



' E N ■ 



(12) 



This Hamiltonian constitutes the starting point for our 
subsequent analysis. It can be solved by diagonalization 
of the (classical) problem of the boson oscillator q = 1 
coupled to the bath oscillators (see Appendix [S] and be- 
low, also Ref. ll5|L Note that coupling of a Luttinger 
liquid to a linear bath has been considered in Ref. 0, 
although the physics discussed there (as well as the cal- 
culation) is quite distinct from our model. 

At the end we are interested in quantities relating to 
the fermions themselves, e.g. the occupation numbers 
(elfin) or t ne Green's functions, like (c„(i)c}j(0)). This 
means we have to go back from the boson operators b q to 
the fermion operators, by employing the relations which 
are also used in the Luttinger liquid. In order to do that, 
we have to introduce auxiliary fermion operators ^(rr): 



ip(x) 



-.e iqx b a 



4> = + <p' 

r = e ~^>® /2 . 



(17) 

(18) 
(19) 



The exponential exp(i(p'(x)) in 1)16(1 may be recognized 
as creating a coherent state with the eigenvalues of b q 
prescribed by Eq. 1)15(1 . The other terms are necessary 
to give the correct normalization and phase-factor, and to 
deal with states other than the A-particle ground state. 

The second equality in i(16|) follows from the Baker- 
Hausdorff identity. The "Klein factor" K is defined to 
commute with the boson operators bq and to produce 
the noninteracting (N ~ l)-particle ground state out of 
the noninteracting A-particlc ground state. Its time- 
evolution follows from [K, H] = [K, E$] = Ku a {N - 1) 
as K(t) — K exp(— ioj (N — l)t). The factor X(x) is given 
by exp(i(A - \)x)/y/5H. Actually, [tf,<p] = - £JL X l/q 
diverges, so we would have to introduce an artificial cut- 
off e~ aq (a — > 0) into the sum. However, this drops out 
in the end result, because r from Eq. I(16|) is canceled 
by the contributions from the equal-time ^-correlators in 
the exponent of Eq. I(21|) . 

Using the relation 1(14[) . we have, for example: 



1 



2ir Jo 
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e~ mx ip(x) dx 



(13) 
(14) 



Note that the ip(x) are not directly related to the fermion 
operators of the particles in the oscillator (which would 
involve the oscillator eigenfunctions) . We have effectively 
mapped our problem to a chiral Luttinger liquid on a 
ring with a coupling oc Fcos(x) (x € [0, 2tt[), see Fig. Q] 
(right). The following results also describe relaxation of 
momentum states in that model. The ip operators are 
useful because they fulfill 



b q ^{x) =-^- e - i ^i>(x 



(15) 



This means the application of ifi(x) on the noninteracting 
A^-particle ground state creates an N — 1-particle state 
that is a coherent state with respect to the boson modes, 
i.e. an eigenstate of b q for every q. As a consequence, the 
fermion operators ip(x) may be expressed asi^: 



ijj(x) = KX(x)e 



itpi (x) e i<p(x) 



KXe^r, 



(16) 



(cl(t)c n ) =—l e^C"'-*) (ft(a/,t)${x,0)) olxdx 1 . 

(20) 

The Green's function involving tp may be expressed di- 
rectly in terms of the correlator of the boson operator <f>, 
using Eq. ifT^jl : 



e -i<j>(x',t) e i<l>(x,0)} r 2 



(21) 



We have used the abbreviation np = N — 1 for the quan- 
tum number of the highest occupied state (in the nonin- 
teracting system). The expectation value on the right- 
hand side is evaluated using the Baker-Haussdorff iden- 
tity and the well-known Gaussian property of the bosonic 
variables (which is unchanged by the coupling to the lin- 
ear bath). This yields: 



exp 



( 22 ) 

These results permit us to calculate the hole propagator 
(cjj(i)c n ), from which we obtain the equilibrium density 
matrix by setting t = 0. Note that we have particle-hole 
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symmetry in our problem, such that the particle propa- 
gator gives no additional information. Writing down the 
expressions analogous to (12011 and l|21|l. and using the 
properties of A(x) and K(t) defined above, we find: 



e--o(^+D* . (23) 

It remains to calculate the correlator of <f> in the interact- 
ing equilibrium: This is where the coupling to the bath 
enters, since the original q — 1 mode will get mixed with 
the bath modes. All the other boson modes are unaf- 
fected. The different boson modes remain independent. 
We obtain (using Eqs. (|T7|> and jTBjl): 



q=l 



[e*^b q {t)-h.c){e^ B b q -h.c)) . (24) 



The expectation values for q > 1 are the original ones, 
for example: 

(&,(<)&*} = e-^(n(wog) + l), (25) 

with n(e) = (e' 3 ' 5 — 1) _1 the Bose distribution function. 

In order to obtain the time-evolution and equilibrium 
expectation values relating to the q = 1 mode (cm. 
mode), we have to diagonalize a quadratic Hamiltonian 
describing the coupling between this mode and the bath 
oscillators, which is described in Appendix lAl 



IV. EVALUATION OF THE GREEN'S 
FUNCTION FOR T — 

At T — 0, some simplifications apply to the evaluation 
of the </>-correlator and, consequently, the hole propagator 

(cUt)c n ). 

In Appendix lAl we express the desired cm. correlator 
in terms of the uncoupled normal mode operators. The 
result can be written in terms of the eigenvalues Clj and 
corresponding eigenvectors C.f. 



6i(t)6] 




(26) 



By defining the spectral weight W(cu) of the cm. mode 
(entry in this notation, see App. 0) in the new eigen- 
basis, 



which is normalized to 1, we can rewrite i|26|l as: 



W(u)(— )e~ lu)t dto = ait) 



In a completely analogous fashion, we find 



&J(t)Si\ =/?_(<) = 



-iut 



(28) 



(29) 



(30) 



1 

4 Jo 



W(uj) 




' du, 



and (b{{t)b{j = (&i(*)6i ) = a(t). 

The function W{uS) is derived in Appendix lAl in terms 
of the bath spectrum. The result is 



W{UJ) = ^ H (^_ a ,2_ A( ^ ) 2 + (r(w 2 )/2) 2 



with: 



A(e) 



2tt 
1 

2^ 



N 




T=0 




(pp) 




m 










du . 


I 




e — v 





(31) 

(32) 
(33) 



(a principal value integral is implied in the last line). 

For weak coupling, W(u>) is a Lorentz peak centered 
around loq. That is why a{t) and P-{t) are small, since 
the spectrum W(lo) is multiplied by a function that van- 
ishes at the resonance near ujq. Both a and /3_ vanish 
exactly at zero coupling. In contrast, /3+(£) describes 
damped oscillations around ujq (see below for a discus- 
sion of the damping rate), starting from /3+(0) 1. 

Evaluation of the cm. mode contribution to the corre- 
lator of 4> (given in Eq. I|24ll ') then leads to the following 
expression: 



y x 6i(t) - ft.c.)(e w 6i -h.c 

(e l(x ' +x) + c.c.)a{t) - e 1 ^' - x) (3+{t) - e i{ -'~^ P- (t) .(34) 

Now we are prepared to evaluate the Green's functions 
of the fermions in the damped oscillator. It is convenient 
to introduce the abbreviations: 



Nb X = e lx , X' = e lx ' . (35) 

^(w)=X;Cg/(w-n i ), (27) 

p~ Then Eq. $Q leads to: 



6 



where 



^(/3 + (t)- e --"*) + ^/3_(t)- 



X 1 

a(t)(XX' + {XX'y 1 ), (36) 



where the subscript (0) refers to the case without cou- 
pling to a bath: 



(0) 



g {e-^x>/xy ^ (3?) 



9=1 



We insert this expression for the ^-correlator into the ex- 
ponent \\'l'2\i that appears in the ip Green's function, Eq. 
B21JI . and perform the Fourier transform (|2t)|l with re- 
spect to x, x' to go back to the original oscillator fermion 
operators c n . 

Then we find that the hole propagator 
i(t)c nF -&n) is given by the prefactor of 



'tif — 8 
~i8n' ! v8 



X n jX in the power series generated from 



J2(X'/X)\ 



,iu> (n F -k)t SE(X ,X' ,t) 



(38) 



fe=0 



with 



SE(X, X\ t) = l- /3+(0) - /3_ (0) + 



«(0) 



X 2 +X '-2 + X~ 2 ) - 

{X'X + 1/(X / X))a(t) + 



(X'/X)(P+(t)-e 



-iiJot\ 



(X/X')0-(t). 



(39) 



Equation i)38|) constitutes our central exact result for the 
single-particle Green's function of the bosonized model. 

The exponent 5E describes the influence of the bath 
on the Green's function. It follows by inserting l|36|) 
into the exponential Eq. (|22|) . The non- interacting 
value of the exponent has been split off and is accounted 
for by the prefactor in Eq. I|38[l . which reproduces 
the correct result for the non-interacting case, namely 



— gi(n f 



-&n)u> t for g n > q_ 



The weak-coupling approximation consists in setting 
a(t) = /?_ (t) = 0, f3 + (0) = 1 and keeping only the (slowly 
decaying) f3+(t), which is approximated by /3+(t) rj 
e -iuj t-Ntt/2 _ jj ere w ^ j s shifted with respect to the fre- 
quency ujq, being equal to woo m the limit ui c 3> u>q (see 
Eqs. H49|) an< i tHDjl f° r the special case of the Ohmic 
bath). In this approximation, we obtain: 



JuJ nt 



rn—0 



ml 



(40) 



v{t) = exp(-i(w - oj )t - N^t/2) - 1 



(41) 



with the decay rate Nj evaluated for the Ohmic bath. 

It is even possible to derive a closed integral expression 
for the exact single-particle Green's function. We only 
state the result: 



1 + (—1)™ n ifn+n'W 



(n+n')w t/2-2/3_(0) 



47T 

i(n F ~ (n' +n)/2)x 

dx ■ 



1 + 0+ - e" ra 

e [e i »«'Wt)-l] e - i '+e- il '» , /S_(t]e il x 

/»/_„ [2a(0) cos(x - uj t) - 2a(t)]. (42) 

Here I m is the modified Bessel function of the first kind. 
A useful approximation becomes possible in the limit 
n <C tif, i.e. for a high excitation energy, even at strong 
coupling: 

exp(e^°*/3 + (i) - /3+(0) + e -' iWot /3_(<) - 0_(O)) x 

e lbJo{nF - Sn)t Io[2a(0)cos(u t) - 2a(t)] . (43) 

Note the trivial dependence on 6n in this limit (a shift in 
frequency). As Iq(x = 0) = 1, we can easily recover the 
weak-coupling form from above. 



V. DISCUSSION OF EQUILIBRIUM DENSITY 
MATRIX AND GREEN'S FUNCTION 

Setting t = in Eq. PS)) yields the equilibrium density 
matrix p nn ', for which only a(0), /3_(0) and /3+(0) = 1 + 
/3_(0) are needed. For the Ohmic bath, and in the limit 
lu c max{iV7, wqoIj we have found (with S = Nj/uqo) 



a(0) 
/M0) 



Li 



W ^00 / ^ 

woo wo \ 2 
w _ woo _ 
woo wq \ 2 



L 2 (44) 
1 



(45) 



where Li = In ^ and 



In 5 2 /2-1 + <5v/<5 2 /4-1 



27rV^2 



(46) 



The behaviour of the density matrix is shown in Fig. 
12 Interestingly, the same form of the equilibrium den- 
sity matrix has been obtained in a model of interacting 
fermions inside a harmonic oscillator-, with a special in- 
teraction, and without bath. In particular, Eq. (43) of 
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100 



n Pun' 




Figure 2: Equilibrium density matrix. Left: Populations p nn 
as a function of coupling strength. Right: Full density matrix 
p nn i for fixed N^y/ujoo = 10, magnitude indicated by radius, 
negative values indicated by different color/gray level. Note: 
Off-diagonal elements have been enlarged by a factor of 3 for 
clarity. We have chosen oj c /woo = 10 4 in both plots. 



Ref. may be compared^ directly to our Eq. I|42|) f° r 
the Green's function, evaluated at t = 0. Naturally the 
dynamics of the two models is different. 

As evident in the plot, there is a 'jump' in the distri- 
bution at n = rip, which is given by: 



bp 



* dx 
2^ 



e -2/3-(0)(l-co S (x)) /o [ 2Q ,( )( cos ( a .) _ jyj 



(47) 

In the strong-coupling limit, the shape of the distribution 
is approximated by a 'discrete error function', 



\ y 



a -fe 2 /(4/3_(0)) 



(48) 



Regarding the Green's function, the series expansion 
of Eq. I)38JI has been carried out using a symbolic com- 
puter algebra system. In order to obtain the Green's 
function in frequency space, the resulting terms of the 
form a(t) na (3 + (t) n + j3-(t) n - have then been Fourier- 
transformed (which leads to repeated convolutions of 
their Fourier spectra, i.e. essentially W(oj), see Eqs. 
P8 ^ .lj59 jl .ffl ty ). This is done numerically (using a fast 
Fourier transform). The results obtained for the Green's 
function in the case of coupling to an Ohmic bath are 
depicted in Figs. 0and0] (In the expansion of the expo- 
nential, the maximum combined power of a, f3± that was 
kept has been 6 for Fig. [3] and 10 for Fig. 

The shape of the broadened peaks in the incoherent 
background reflects the Fourier-transforms of a(t),/3± (t). 
All of them are related to the spectral function W(oj) de- 
scribing the damped cm. motion (see Eqs. 1281 1291 and 
I5U) |. In particular, at weak coupling, the Fourier trans- 
form of the dominant contribution (3+(t) is essentially 
given by W(oj), see Eq. (j55J). 

For the Ohmic bath, whose power spectrum is given in 
Eq. ©, we have (see ©): 



-Do 



+ -Njoj c 



(49) 



The resulting spectral function W(uj) associated to the 
damped cm. motion, sec Eqs. (|2*7|) and HA17jl . is found 
to be 



W(lu) = ^oj 2 x 

7T 



(50) 



2 2 jN 0J c + UJ 

UJ — Wnn UJ ill 



7T \oj c - U)\ 



(jNuj) 2 



This expression holds for < uj < uj c . (The denominator 
becomes at some uj > uj c , which yields a 8 peak whose 
weight vanishes for 7iV — > 0. In order to obtain this 
peak, one has to replace 2jNu> by an infinitesimal +0 for 
lj > lo c . This would be different for smooth cutoffs, and 
it does not affect any of our results). The function W(ui) 
is displayed in the inset of Fig. [3] (where uj c /ujqo = 10). 
At low frequencies uj -C uj c and damping rate 7^V <C 
u) c , W(uj) is proportional to the magnitude squared of 
the susceptibility derived for the linear response of the 
velocity of the damped cm. oscillator, with renormalized 
frequency luqq and damping rate jN: 



\XvA")\ 2 « uj 2 [{uj 2 - uj 2 f + ( 7 jVw) 2 



(51) 



In the overdamped case, > 2woo, the zeroes of the de- 
nominator occur at purely imaginary frequencies. How- 
ever, the shape of the function (|51|l shows a qualitative 
change already at jN — y/2ujoo, when the curvature at 
uj = changes sign (in the susceptibility the maximum 
vanishes). This may be observed by comparing the low- 
est curve in Fig. 0] corresponding to ^N/ujoo = 1.6, to 
the other curves. 

Apart from the generic exponential-type decay at in- 
termediate times, we can also evaluate the behaviour of 
a(t), P+(t) and /?_(£) at long times t » l/min(u;oo, ^7)> 
for the special case of the Ohmic bath. All three func- 
tions then decay like t~ 2 , namely as a(t) sw (3±(t) « 
— [Njujo/C^ttujqq)]/^ , which corresponds to the decay of 
the incoherent part of the Green's function. 

The general long-time behaviour of the Green's func- 
tion (independently of the bath spectrum) can be read 
off directly from the exact expression given above in Eq. 
P%)l: Since the coupling to the bath damps away the 
cm. motion, the functions a(t), (3+(t) and /3-(t) all de- 
cay to zero at t — > oo (excluding cases such as a gapped 
bath) . The exponent SE in Eq. i)39|) thus approaches a fi- 
nite value. This leads to the conclusion that the Green's 
function (cl lF _ Sn (t)c np -s n ^ does not decay to zero in 
the limit t — > oo, for any 5n > 1, since contributions 
from (X' /X) Sn with Sn > 1 remain. In frequency space, 
it consists of a number of delta peaks superimposed onto 
an incoherent background related to the damped cm. 
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Figure 3: Green's function of fermions in a damped harmonic 
oscillator (with an Ohmic bath): Fourier transform G nn (oj) 
of (cjj(i)cn) plotted vs. uj, for different values of n (curves 
displaced vertically for clarity) . Note that u> is measured with 
respect to — tifojq and we have set o>o = 1. The strength of the 
coupling is given by N^f/ujoo = 0.4. Inset: Weight function 
W(u) vs. a). The height of the smaller S peaks is an indication 
of their weight. 
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Figure 4: Dependence of Green's function on coupling 
strength: Fourier transform G„„(ui) of (c^(t)c„) vs. uo, for 
different values of the coupling strength, indicated in the plot 
(curves displaced vertically for clarity). Note that u) is mea- 
sured with respect to —uf^o and plotted in units of luq, and 
that too changes with 7, for fixed u>oo, see Eq. 149H . This plot 
is produced for a fixed excitation number, n — uf — 4. 



mode (see plots). These facts can also be seen from the 
simple weak-coupling expression, Eq. (I4UII . 

This result may be unexpected, since according to a 
naive application of the Golden Rule to the Hamiltonian 
|T)l. any electron inserted into the system at n > tif + 1 
(and any hole inserted at n < np) should decay towards 
lower (higher) single-particle energy levels by sponta- 
neous emission of energy into the bath (for T = 0). The 
rate is expected to be: 



GR 



2 Jpp) J*. 



(52) 



(with n ~ riF ~ N). For an electron inserted just above 
the Fermi surface (n = uf + 1), or a hole inserted at 
n = np, the Pauli principle blocks relaxation and the 
corresponding Golden Rule rate vanishes. The decay 
rate of Green's functions is expected to be Tqx/2. In 
the following, we discuss the deviations from this naive 
expectation. 

For weak coupling, the decay of the Green's function 
(cjj(i)c„) calculated above is determined by the decay 
of P+(t), while a and /3_ are approximately zero. For 
example, at n — np — 1 the hole-propagator is directly 
proportional to P+(t): 



iriFUjQt 



(53) 



The decay of this function is determined by the width 
of the Lorentz peak appearing at u) w luq in the "spec- 
tral weight" W{ui) defined above (see Eqs. (|27Jl and 
(I29|l ) and calculated in Appendix 1X1 According to Eqs. 
(IA17|I and (|A15|1 . at weak coupling this width is given 
by Auj — T GR /2, as expected. Therefore, there is quan- 
titative agreement between the exact calculation and the 
Golden Rule result at n = np — 1. This is true also 
at n — rip, where the Green's function evaluated in the 
preceding section only has a delta peak, related to Pauli 
blocking (apart from corrections that vanish in the limit 
of weak coupling), see Fig. |3| However, at n < np — 1, 
there is always an additional S peak that is incompatible 
with the Golden Rule expectation, as explained above. 
Essentially, this discrepancy stems from the fact that 
the many-particle states in the harmonic oscillator are 
in general strongly degenerate. Therefore, the precondi- 
tions for the usual derivation of a master equationil are 
violated. We will trace back this behaviour to the fact 
that the exact wave function splits into cm. and rel- 
ative motion (see the remainder of this section and the 
following section) . Readers not interested in this detailed 
analysis may wish to skip directly to Section lYlII for the 
time-evolution of the density matrix. 

The excited many-particle state | , F), which is created 
by adding an extra hole (or extra particle) and whose 
time-evolution determines the Green's function, can be 
written as a superposition of states, each of which dis- 
tributes the given number of excitation quanta Sn differ- 
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ently onto center-of-mass ("cm") and relative ("r") mo- 
tion. Therefore, in the weak-coupling limit, where one 
starts out of the factorized ground state |0) cm |0) r |0) B , 
it is of the form: 

|*> = |0) B (oo |0) cm \6n) r + a, |l) cm \Sn -l) r + ...). 

The first state does not contain any extra excitation of 
the cm. mode. Therefore, it is unaffected by the c.m.- 
bath coupling and does not decay. This leads to the 
main 5-peak contribution at a frequency Snuo, which is 
determined by the number of excitation quanta Sn. 

The other states will decay. The resulting broadened 
peaks are shifted towards lower frequencies, due to the 
renormalization of the cm. mode frequency to the value 
ujoq. Thus these peaks occur at frequencies of the form 
(Sn — n cm )ujQ + n cm u>oo, where n cm counts the number of 
quanta in the cm. mode. This is visible in Figs. [21 and 

where the Fourier transform of the Green's function is 
displayed. In the limit of weak coupling (a — (3- — and 
/3+(0) = 1 ), the expansion yields the following weights 
of the different peaks: 

! *-=sr ( _i)i 

"■cm- j=0 J- 

for < n cm < Sn. These weights are normalized 
(Sn™ =o l a »c m | 2 — 1); since the application of a parti- 
cle annihilation operator onto the non-interacting ground 
state creates a normalized state |^) (for Sn > 0). In par- 
ticular, the weight of the n cm = 0-component, which 
leads to the S peak, goes to the constant value of 1/e 
in the limit Sn — > oo. Thus, there is a non-decaying 
contribution even at arbitrarily high excitation energies. 
This fact is connected to the assumption of large N un- 
der which the present results have been obtained (see the 
following section). The weight of the contribution from 
n cm = Sn — 1 always vanishes, as may be observed in the 
plots as well. 

For stronger coupling (Fig. 0}, other S peaks start to 
grow both above and below the frequency of the main 
peak. These arise because the interacting ground state 
contains contributions from excited cm. states and bath 
states. We can write it down in the following schematic 
form, where \j) B denotes a bath state for which the sum 
of all harmonic oscillator excitations is equal to j: 

\GS) cm+B = (56) 
coo |Q) cm |0) B + c n |l) cm |l) fl + c 20 |2) cm |0) B + . . . . 

Any of the excited states discussed above, containing n cm 
cm. excitations (see Eq. 1154(1 ). may therefore have a fi- 
nite overlap with this ground state of the coupled system 
cm. -bath. Thus, a non-decaying contribution at the fre- 
quency (Sn — n cm )oJo of the remaining excitation in the 



relative motion survives. This explains the S peaks at 
frequencies equal to integer multiples of ljq, situated be- 
low the unperturbed excitation energy of SnajQ. There 
are also peaks at still larger frequencies, because the ad- 
ditional particle (or hole) is introduced into the many- 
particle ground state that already contains excitations of 
the cm. mode, see Eq. I|56f) . Therefore, the resulting ex- 
cited many-particle state |W) will also contain contribu- 
tions with more than Sn excitations of cm. and relative 
motion combined, in contrast to the weak-coupling form, 
Eq. (|54|l . Since adding Sn excitations to one of the states 
\j) cm |0) r that appear in the many-particle ground state 
may lead to states containing up to Sn + j excitations 
in the relative motion alone, this explains the appear- 
ance of S peaks at (Sn + j)u!o in the Green's function. 
For the same reason, the "incoherent background" due 
to the cm. mode extends beyond the main frequency of 
SnajQ. 

Finally, it may be noted from the figures that the ad- 
ditional 8 peaks appear only at frequencies of the form 
(Sn + 2j)u>o, i.e. they are removed by an even multiple 
of luq from the main peak. Physically, the reason is the 
following: The coupled ground state of cm. and bath 
remains symmetric under inversion of the cm. and oscil- 
lator coordinates. Therefore, the sum of cm. and bath 
excitations is always even, as indicated in (|56|l . When 
adding Sn extra excitations to the many-particle ground 
state (which is (|56H . multiplied by the ground state of 
the relative motion), a part of those excitations may go 
into the cm. motion, leading to peaks below and above 
the main frequency. Only if the number of added or sub- 
tracted cm. excitations is even, a nonvanishing overlap 
with the many-particle ground state may develop, lead- 
ing to a non-decaying component in the Green's function. 

The considerations of the current section apply to the 
approximate Hamiltonian l|12|l that is good in the limit 
N — > oo of large particle number, and which has been 
amenable to bosonization. However, the qualitative ar- 
guments concerning the contributions of cm. and rela- 
tive motion remain valid for finite N as well. This is the 
topic of section IVll 



VI. FINITE PARTICLE NUMBER: 
CENTER-OF-MASS MOTION IN EXCITED 
FOCK STATES 

In this section we analyze in more detail the splitting 
into cm. and relative motion for arbitrary finite N. We 
will explain that the weight of the coherent 5 peak in 
the Green's function decays when moving towards higher 
excitation energies, on a scale set by the number N of 
particles. This is why, in the limit N — > oo considered 
for the bosonized model, this weight even saturates at 
a constant value for excited states arbitrarily far above 
the Fermi surface. We will not actually evaluate Green's 
functions for the finite N case, which is still much more 
difficult. 
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For brevity, we set m — loq — 1 throughout this section. 
The coordinates and momenta of the individual particles 
are related to the harmonic oscillator raising and lowering 
operators by 



1 ,„ , U 

T2 {aj + aj) 



—;=(aj 
V2 J 



a}) 



(57) 



where it is understood that a,j acts only on the coordinate 
of particle j. The center-of-mass motion is described by 



1 N 

3=1 



N 

E^ 

3=1 



(58) 



such that we obtain the following operator that lowers 
the excitation of the cm. harmonic oscillator by one: 



A 



N » P 
2 {X + 'N ] 



N 

E« 

3=1 



(59) 



N. 



Here we have used that the total mass is M = Nm 

Let us consider an iV-particle Fock state \^N.8n) con- 
sisting of the Fermi sea filled up to (and including) level 
rip — 1 = N — 2 and an additional single particle that has 
been placed into the excited level at tif + Sn. It is our 
goal to find the probability P n of having n excitations in 
the cm. mode, given this many-particle state. In partic- 
ular, the probability of having excitations will be the 
weight of the "coherent component" that does not decay 
in spite of coupling to the bath. This holds quantita- 
tively in the weak-coupling limit, where we may neglect 
the change in the eigenstates of the cm. mode, see the 
previous discussion. For arbitrary coupling, it is still the 
correct qualitative picture. 

First, we rewrite A in second quantization: 




Figure 5: Weight Po of the center-of-mass ground state (Eq. 
1651 '! as a function of excitation Sn, for different particle 
numbers N, including the bosonization result from Eq. (15511 
("N = oo"). 



Now we make use of the fact that the total number of 
excitations in the cm. mode is limited by Sn, i.e. P n = 
for n > Sn. Therefore, it is possible to start at j = 
Sn, equate (|61fl to (|62[1 . solve for Ps n , and then proceed 
iteratively all the way down to Po- In each step, only 
the probabilities that have been calculated before appear 
in addition to the unknown Pj. We can write down the 
set of equations very transparently by introducing the 
abbreviation 



^ = ^- l ^ WT^ NSnp ^- (63) 

Then we have (from equating l|61l) to H62H - ) the recursive 
relation: 



oo 



(60) 



Using this, it is straightforward to check that 



np 



N 



5n + np — j + 1 

N ' 

(61) 

for j < Sn, otherwise this is 0. On the other hand, if we 
think of \^N,5n) as being written in a basis that splits 
relative and cm. motion (compare Eq. §54$ ) and take 
into account the usual matrix elements of a harmonic 
oscillator lowering operator, we obtain for the same ex- 
pectation value (with P n the occupation probabilities for 
the different cm. states): 



oo 

E 

71= j 



P n n(n - 1) 



■i + i). 



(62) 



N 1 ' 



Pi 



(N + l-1) 



(N + l) 



E 

fc=i 



fc! 



(64) 



The solution starts with po = 1. Note that we always get 
Pi = (for any N), corresponding to Ps n -i — 0: This has 
been observed already in the context of the bosonization 
solution (see previous section) . It is also remarkable that 
the equation for pi does not depend on Sn. Therefore, 
we can solve the problem at once for arbitrary Sn (which 
then only enters in P n via Eq. (|63|I % 1. After solving these 
equations, it is found that the weight of the cm. ground 
state, 



decays as a function of excitation energy Sn (for large 
Sn), on a scale set by N, see Fig. |SJ 
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Finally, we remark that the general structure of the ex- 
act eigenf unctions of our problem (defined by the original 
Hamiltonian, Eq. 0J) remains identical to that found in 
the limit N — > oo, using bosonization. For arbitrary 
system-bath coupling, the eigenstates of the full system 
will be of the form 



CM-B 



) ® \V el ) , 



(66) 



where "SB" refers to the full Hilbert space of fermions 
coupled to a bath, while \<fr CM ~ B ^ is a new eigenstate 
of the coupled system cm. /bath and \ip rel } is an unal- 
tered, antisymmetric eigenstate of the relative motion. 
However, the evaluation of quantities like Green's func- 
tions or the reduced density matrix, starting from this 
expression, becomes very cumbersome, which is the rea- 
son why we have restricted most of the discussion to the 
bosonization solution. 



VII. TWO-PARTICLE GREEN'S FUNCTION: 
DECAY OF POPULATIONS AND DEPHASING 

While the single-particle Green's function reveals how 
fast an electron is scattered out of an initial state, we 
have to turn to the two-particle Green's function in order 
to learn about the time-evolution of the reduced single 
particle density matrix. This describes relaxation and 
dephasing. For our purposes, we will be particularly in- 
terested in the following two-particle Green's function of 
fermions in the damped oscillator (the evaluation of other 
four-point correlators proceeds analogously): 



c n (0)clm(t)cl(0) 



(67) 



Without coupling to the bath, we could apply Wick's 
theorem to obtain: 



;e„(o)4(t)Q(t)&(o); 

Sni'S w (l - (n„))(l - (rv»e iw °("-"')' + 

5nn'5iv(l-{n n )){hi) , (68) 

If the bath damps the motion of the fermions, we may 
rewrite (16711 in terms of the auxiliary fermion operators 
ip(x) defined in (|14H : 



1 



(2^) 2 Jo 



dxdx'dydy' e *C*'+nV-i*-n») x 



il>(y,0)ft(x',t)i>(x,tW(y',0)} (69) 



Starting from the representation of ip in terms of boson 
operators <p, Eq. (|16|) . this may be evaluated using the 



same methods as in Sections IIIII and IIVI In Appendix 
IB1 we show how the two-particle Green's function is ob- 
tained via a series expansion, similar to the single-particle 
Green's function. The general exact result is given in 
Eqs. ljB7(l - l|B9() . Although in principle this expression 
permits evaluation of the two-particle Green's function 
at arbitrary coupling strength, this task has proven to 
be computationally much more difficult than the analo- 
gous calculation for the single-particle Green's function, 
due to the larger number of terms generated in the power 
series. Nevertheless, the weak-coupling results already 
contain many nontrivial features. 

We may now answer questions like the following: If we 
introduce an extra single particle in some level n above 
the Fermi sea, how does it decay towards lower-lying lev- 
els, by emitting some of its energy into the bath? We 
first look at the time-evolution of the populations 



Pnn{t) = (c S 4(*)Cn(i)4 



(70) 



where the two outermost operators create the desired 
state at time 0, while the two operators in the middle 
test for the population. We note that, in principle, the 
newly created state may have norm less than 1, if the 
level n happens to be partially occupied already in the 
ground state on which the creation operator c~ acts. In 
that case, it would be necessary to divide iffO)) by the 

norm, ^CfictV However, since we will restrict the evalu- 
ation to the weak-coupling case, the states above np are 
initially empty (at T = 0). 

In particular, we have been able to derive simplified 
expressions in the limit of weak-coupling and high initial 
excitation energy (see Appendix 0. We obtain: 



^ ri (0)4(t)Q(t)<t(0)) « e^C™-™')*^,^^, x (7i) 

{/Odccay(™' ~l,t) + pheatQ ~ Tip - 1, 1' - Up - 1, t)} 

with the first part describing the decay of the excitation, 



Pdccay (m,t) 



(-iy 



{v{t)+v*{t)) m e v ® +v *® , (72) 



where m = hq — n = n' — n' may be interpreted as the net 
number of quanta transferred to the bath (v(t) is defined 
in Eq. lUTl) ). The other part describes "heating" around 
the Fermi surface (see below): 



Ph ea .t{n,n',t) = i/(t) r 



2(mi+m 2 ) 



^-^ mi!m2!mi!?7i2! 



(-1) 



m.2+rn2 



(73) 

where the triple sum runs over rh\ = max(0, n' + l) . . . oo, 

7112 = . . . OO, TTll = max(0, 7772 + 77 + 1) . . . 77- — 7i' + 777l + 

7772i and we have 7772 = 777! + 7772 — Til +77 — 77'. We note 
that the limiting case (|71|l appears to be a very good 
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Figure 6: Time-evolution of phoat (n,n' , t), describing heating 
around the Fermi surface due to an extra particle placed in 
a highly excited state above the Fermi sea (Eq. |JZ3)- Sec 
text for connection to two-particle Green's function, Eq. I71t . 
In this plot ujq — uj' . The radius of each dot gives |p n „/(t)| 
(being 1 for the two occupied states in the lower left corner). 
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Figure 7: Time-evolution of the populations p nn (t) (Eq. 1701 1 
for an electron placed in an initial state h = np + 6 above 
the Fermi sea (weak-coupling result) with u>o = ui' (left) or 
ljo = u>' +2Ny (right), see text. The initial and final states are 
the same in both cases. Different curves have been displaced 
vertically for clarity. Note the incomplete decay and heating 
around the Fermi surface. 



approximation to the full result, even for small excitation 
energies. The behaviour of p nn (t) shown in Fig. gives 
an impression of /?decay(w, i) for the levels not directly 
near the Fermi surface (with m being the distance from 
the initial excited state). A plot of pdccay can be found in 
Ref. [H In Fig. H 

we have shown the time-evolution of 
Phcat, which results in a good approximation to the full 
evolution displayed in Fig. [SJ taking into account Eq. 

CD- 

The time-evolution of the populations, p nn (t), is shown 
in Fig. which has been calculated in the weak-coupling 
limit, but without the approximation of a high initial ex- 
citation. Note that, in general, the difference between u>q 
and us' can be of any sign, depending on the details of 
the bath spectrum. This frequency shift between the cm. 
mode and all the other boson modes may lead to beating 
oscillations that may be visible in the time-evolution of 
the populations or the density matrix, which is shown in 
Fig. [3 right. In Figs. (left) and [SI however, we have 
assumed lu'q — loq (e.g. the bath spectrum is symmetric 
around the transition frequency, for a fixed decay rate of 
-/V7). We note that, contrary to naive expectation, the 
particle does not decay all the way down to the lowest 
unoccupied state np + 1. Rather, in the long-time limit 
(which is already reached at A^i w 5 to a good approx- 
imation), the extra particle is distributed over the range 
of excited levels above the Fermi surface, up to the ini- 
tial level h. Again, this is because only the cm. mode 
couples to the bath, such that, even at T = 0, a fraction 
of the initial excitation energy remains in the system. 
For the limiting case of high excitation energy, discussed 
above, this may be seen from pdecay^j t) — > (2 m /m!)e~ 2 . 
Moreover, we see that the population of the highest oc- 
cupied states in the Fermi sea decreases. The fermions in 
these states become partly excited by absorbing energy 
from the extra particle, due to the effective interaction 
mediated by the bath. This is described by pheat of Eq. 

Another application of the two-particle Green's func- 
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Figure 8: Time-evolution of the single-particle density ma- 
trix p nn '(t) (Eq. 1741 1 starting from an extra particle in a 
superposition of states ni = 5, 712 = 6 above the Fermi sea 
(u)q = Wq). Levels n,n — uf — 1, np, iif + 6 are indicated 
by grid lines. 



tion (|B1|) consists in looking at the decay of a superpo- 
sition of states. We imagine that, at time 0, a particle is 
placed in an equally- weighted superposition of two levels 
n\ and 712 above the Fermi surface. This will introduce 
off-diagonal contributions in the single-particle density 
matrix of the fermions. The subsequent relaxation will 
suppress the coherence (i.e. the off-diagonal elements) 
and transfer population to the lower levels. In Fig. |SJ 
we have plotted the time-evolution of the density matrix 
defined by 



'(*) 



(c ni + c*,)4(*)S»(*)(4i + 4 2 )) . (74) 
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Again, we have evaluated only the weak-coupling case, 
such that no additional normalization factor is necessary 
(for ni, n% > rip). 

Apart from the features already mentioned above, we 
observe the decay of the off-diagonal elements (i.e. de- 
phasing) to be incomplete. A part of the coherence sur- 
vives in the relative motion which is unaffected by the 
bath. This is also evident from the limit pdccay {m, t — > 
oo ) discussed above. If a slight anharmonicity were in- 
troduced in the potential, such that cm. and relative 
motion are no longer independent, we would expect to 
see the same behaviour at short to intermediate times. 
Only in the long-time limit, the fermion system would 
relax fully to the N + 1-particle Fermi sea ground state. 



VIII. RELAXATION WITH ARBITRARY 
ENERGY TRANSFER 

Up to now, we have considered a coupling be- 
tween system and bath where the particle coordi- 
nates enter linearly (since this corresponds to a direct 
many-particle generalization of the well-known Caldeira- 
Leggett model). The coordinate operator only couples 
adjacent oscillator levels. Therefore, a particle intro- 
duced somewhere above the Fermi sea will have to relax 
in a series of spontaneous emissions, each of them leading 
down exactly one level. For this reason, the influence of 
Pauli blocking becomes apparent only immediately above 
the Fermi level. In addition, this type of coupling also 
leads to damping of the center-of-mass motion only, such 
that the system is non-ergodic. 

In this section, we are going to discuss a more general 
model, where the bath may induce transitions between 
oscillator levels that are farther apart. The most impor- 
tant consequence will be that the presence of the Fermi 
sea leads to a strong dependence of decay rates on the dis- 
tance to the Fermi level. As one moves towards the Fermi 
level, more and more of the transitions will be blocked 
due to the Pauli principle. It is interesting to observe the 
consequences of this in an exactly solvable model, where 
we are able to go beyond a Golden Rule description. Fur- 
thermore, this coupling will (in general) lead to ergodic 
behaviour, i.e. full relaxation of the excitation that has 
been added to the system. 

We start from the bosonized Hamiltonian, Eq. (|12|1 . 
and replace the coupling term by a more general form: 

I TV °° 

q—l n 

I — AT 00 

V2^£/.^+*S)- (75) 

q=l 

Now the bath may induce transitions between levels 
n + q and n, with an arbitrary amplitude oc f q (which, 
however, must not depend on n). For fi = 1, f q — 



(q > 1) we would recover the original model. In order 
to simplify the following expressions somewhat, we have 
assumed f q to be real-valued (nothing essential would 
change for complex amplitudes). For f q = const, the 
amplitude of transitions between two levels does not de- 
pend on their distance. The overall prefactor has been 
kept the same only in order to retain similarity with the 
previous results (it may always be compensated by ap- 
propriate choice of the amplitudes f q ). Note that, in 
terms of the original model (before bosonization) , the in- 
teraction given here only approximately corresponds to 
a nonlinear coupling that adds higher powers of the par- 
ticle coordinates, since such a coupling would lead to ad- 
ditional terms in the bosonized version. 

We will take the interaction term Eq. (|75|l as the start- 
ing point of our analysis, which employs the methods in- 
troduced before. In the weak-coupling limit, we expect 
each of the boson modes to be damped with a Golden 
Rule decay rate that follows directly from Eq. J73J): 

r 2 N 2 /ff\ t= ° 

q 2mUjQ q \ I qu) 

Adding a particle above the Fermi sea, in state uf + Sn + 
1, will introduce bosonic excitations of up to q = Sn (in 
the weak-coupling limit, starting from the unperturbed 
system), and its decay will be governed by corresponding 
rates T q . As Sn grows, the typical decay rate will grow 
as well, as expected. This will be analyzed in more detail 
further below. 

Now all the boson modes couple to all the bath oscil- 
lators. We retain the definition of oscillator annihilation 
operators dj given in Eq. ijAlfl for the bath modes j > 1, 
while the boson modes are now defined to have negative 
indices, j = — q < — 1, with d- q = b q and Q- q = quo for 
q > 1 (there is no j — mode in this notation) . Using Eq. 
(|A1|) . this definition leads to y/q{b q + b q ) = qy / 2mui Q- q . 
The resulting problem of coupled oscillators may again 
be written in the form of Eq. I|A2(1 . However, now the 
"perturbation" matrix V that couples boson modes and 
bath oscillators has the following nonvanishing entries 
(j,<?>l): 

After diagonalizing the problem with an orthogonal ma- 
trix C (compare Appendix El Eqs. JA3| l-l |A6|| l). the re- 
quired correlators of boson modes, of the type (b q i (t)b q \ 

can be calculated exactly as before. Now the interaction 
with the bath also mixes different boson modes. There- 
fore, we may have q' =/= q in general. The appropriate 
generalization of the definition for the spectral weight W 
(Eq. (j23) is 

W q >M = J2 C-^jC-^Siu - %) . (78) 



14 



At T = 0, we obtain, using Eq. i|A7|l (in analogy to (J2SJ) 
or (P|> ): 




^+ _ \ du) . (79) 



Analogous expressions hold for the other correlators: The 
sign in the first bracket on the r.h.s., involving q' , is posi- 
tive (negative) if the corresponding operator on the l.h.s. 
annihilates (creates) a particle, the rule for the sign in the 
second bracket is the opposite (compare previous Eqs. 
PS)l. 1(531) . and (|30I) ) - The actual calculation of W q - q {oj) 
with the help of the resolvent is described in Appendix 

The general expressions for the Green's functions still 
look like before (see Eqs. {SJJ, iJHJ, (J2U, (EH) for the 
single-particle Green's function). However, now the eval- 
uation of the correlator (|24f> of will yield contributions 
for all combinations of q, q' : 



= -J2 4p((^ V VW-ft.c)(e l « I i,-/ i .c.)) (80) 



q,q' = l 

Together with the results for the ^-correlators, Eq. (|8*TJl . 
this can be used to evaluate the Green's function, which 
would work exactly as before in principle (expanding the 
exponential in terms of exp(iir), exp(ix')). However, now 
one would have to deal with far more terms. We have not 
carried through this calculation in the general case so far. 

Nevertheless, interesting behaviour of the fcrmions is 
already found in the weak-coupling limit. In this limit, 
we neglect the effective coupling between boson modes 
that has been induced by the bath (see Appendix©! and 
describe the correlator of each boson mode separately as 
a damped oscillation, with: 



b q (0 q 



-i(qujQ-\~Re 8Q, — q )t~Im 8Q- q \ t\ 



(81) 



Here the complex- valued frequency shift Sfl - q is given by 
Eq. 1D9(1 of Appendix IdI In particular, the decay rate is 



I m S^ q =qf^UFF) 



mujQ 2 \ 



qu 



(82) 



As expected, this is equal to one-half the decay rate T q 
given in Eq. I|76|) > which describes decay of populations 
(rather than Green's functions). 

For the special case of the Ohmic bath spectrum, Eq. 
©, (and with constant f q = 1), we obtain a decay rate 
rising quadratically with excitation g, 



iV 7 



Im5£L- q = q } 



(83) 



while the frequency shift is linear in q, up to corrections 
on the order of (qajo/uj c ) 2 : 



Re Sfl- q = 



qNj 

TTUJq 



qio 1 + (quj )/u c 
—u c H In - — — 

2 1 - (qu )/uj c 



(84) 

(Jo 7T 

These expressions hold for quo in the "scaling" region, 
where the bath spectrum rises linearly with frequency 
(in our case, for the sharp cutoff, qcoo < w c)- The last 
approximation requires quio <C w c - 

The linear rise of the frequency shift with q will 
prove convenient in the following calculation, since it im- 
plies that the original frequency uj is just shifted by 
— (LO c /uJo)N"f/ir, for all boson modes with sufficiently 
small q. Concerning the decay rate, the quadratic in- 
crease is expected, since in a naive Golden Rule picture 
the decay of a particle from state np + Sn + 1 should be 
produced by all the transitions from q = 1 up to q — Sn. 

Adding up their rates (which grow like ( FF > oc q in 

\ / q u) 

the case of the Ohmic bath) leads to a total rate oc Sn 2 , 
which corresponds to the decay rate of the highest boson 
mode that is excited by adding this particle. 

We will now confirm this qualitative picture by analyz- 
ing the weak-coupling form of the single-particle Green's 
function for the Ohmic bath, at T — (with f q = 1 up to 
some cutoff). We must assume the weak-coupling result 
to be a good approximation for all q. The approximations 
introduced above yield: 



_ e n(^' -x) e -i^' qt e ~N lq 2 \t\/2 



q =l 



(85) 

Here the shifted frequency is given by lo' = ujq — 
(ui c /uj )N~f /it. Note that, according to the discussion 
above, the simple form of the correlator used here will be 
valid only for sufficiently small q, where the conditions 
qajQ <C u) c , f q = 1 hold. In spite of this, we have ex- 
tended the sum in l|85fl over all q. We have been able to 
do so because only g-values up to the excitation level Sn 
will enter the end-result to be derived below. The condi- 
tion about q being small is therefore a condition on the 
excitation considered. 

By inserting (|85|l into the exponential l|22f) and ex- 
panding it in the usual way, we obtain for the particle 
propagator in this approximation: 

(c nF+1+Sn (t)cl F+1+Sn ) = e^o(n F+ i)t x 

e -iu' Snt y 1 y 1 c -Nf(q2 + ...+ q l)\t\/2oQ} 

^ m! ^ q 1 ■ . . . ■ q m 

m=l gi....<? m 
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Figure 9: Distribution of decay rates Tgf entering the Green's 
function in the weak-coupling approximation given by Eq. 
1861 1 . for the Ohmic bath at T — 0, where the bath connects 
arbitrary oscillator levels. Decay rates are plotted as a func- 
tion of excitation level Sn, and the grey level visualizes the 
weight in Eq. (1861 . The full lines represent lower and up- 
per bounds N'ySn/2 and NySn 2 /2, while the dashed line is 
N-ySn(Sn + l)/4. 



Figure 10: Decay of populations for a bath that induces tran- 
sitions between levels that are arbitrarily far apart (indicated 
by arrows); calculated in the weak-coupling approximation 
explained in the text, for the Ohmic bath at T — (see Sec- 
tion and Appendix lEl. 



IX. CONCLUSIONS 



Here the sum over qi,...,q m is restricted by Ylj=i 1j = 
6n and qj > 1 . Therefore, we obtain contributions from 
decay rates between (N r y/2)5n and (N'y /2)Sn 2 . The 
naive calculation (applying Golden Rule to Eq. (|75|l and 
adding up all the transitions that are not blocked by the 
Pauli principle) would give the following decay rate for 
the population of level np + 1 + 5n: 

N J 2 -. / - ~\ T=0 
2tt- >(FF) = N 1 5n(5n+ l)/2. (87) 

In order to compare with the decay rate of the Green's 
function, we divide by 2 and obtain (Nj/2)5n(5n + 1)/2, 
which is between the lower and upper bounds derived in 
Eq. H86JI . The distribution of decay rates, involving the 
proper weights from (|86Jl . is plotted in Fig. E3 

The same approximation can be used to obtain a 
weak-coupling version of the two-particle Green's func- 
tion (and, thus, the time-evolution of the reduced single- 
particle density matrix, compare Sec. IVII|I . The detailed 
steps of the evaluation are explained in Appendix El An 
example of the resulting time-evolution of the popula- 
tions is shown in Fig. EH The relaxation towards the 
N + 1-particle ground state is complete, the system is 
crgodic. At intermediate times, heating of the Fermi sur- 
face is also observed, although it is much less pronounced 
than in the example of Fig. [7| 



We have analyzed a model system of fermions in a har- 
monic oscillator, that are subject to a fluctuating force 
deriving from a linear bath of oscillators. The force cou- 
ples to the center-of-mass motion of the fermions. This 
is a generalization of the Caldeira-Leggett model of a 
single particle in a damped harmonic oscillator, and it 
is, in principle, exactly solvable. However, here we have 
analyzed the limit of large particle numbers, using the 
method of bosonization, i.e. we have considered particle- 
hole excitations around the Fermi surface. One of the 
boson modes corresponds to the center-of-mass motion, 
experiencing dephasing and damping, while the others 
belong to the relative motion which forms a kind of 
"decoherence-free subspace" . We have derived exact (for 
N — > oo) analytic expressions for Green's functions, that 
can be obtained from the coefficients of a power series 
expansion. These have been evaluated with the help of a 
symbolic computer algebra program. The single-particle 
Green's function has been discussed in detail, showing the 
smearing of the Fermi surface, the effect of Pauli block- 
ing, the dependence of level shapes on the bath spectrum 
and the coupling strength, as well as the appearance of 
coherent peaks, due to the relative motion. Based on 
the expression for the two-particle Green's function, we 
have analyzed the decay of an excited state created by 
adding one particle above the Fermi surface, where one 
can observe the "heating" around the Fermi surface (due 
to the effective interaction between particles), as well as 
the incomplete decay of the excited particle. We have 



1G 



provided simplified expressions for the limiting case of 
high initial excitation energy. In addition, we have dis- 
cussed the time-evolution of the single-particle density 
matrix for an extra particle that has been added in a co- 
herent superposition of excited states. Finally we have 
extended our analysis to a more general coupling, where 
the bath induces transitions between oscillator levels that 
are not adjacent. In that case, the influence of the Pauli 
principle becomes even more pronounced. 



Acknowledgments 

We thank C. Bruder, W. Belzig, H. Grabert, D. Loss 
and A. Zaikin for comments and discussions. The work 
of F.M. has been supported by the Swiss NSF and a DFG 
grant. 



Appendix A: COUPLED OSCILLATOR PROBLEM 

We define the annihilation operators of the bath modes 
(for j > 1) as 



CP 



-(Qj 



ifl 



■) => Qj = (d 



, Tdt)/V2mpy, 
(Al) 

and do = b\, which makes Qo and Pq the canonical vari- 
ables of the cm. mode and the coupling in Eq. i|12|) 
equal to ^/NFQq. 

We now have to diagonalize the following quadratic 
Hamiltonian, according to Eqs. I|12l) and J5J: 



H' = ^-P l P 
2m 



-Q t (Q 2 + V)Q. 



(A2) 



Here P and Q are (column) vectors, with operator en- 
tries Pj, Qj (j = O...Nb), where j > refers to 
the bath oscillators, and j — refers to the q = 1 
mode. The diagonal matrix fl 2 contains the eigenval- 
ues Vt 2 of the original problem (with f2 = uj ) and V 
is the coupling matrix, with the nonvanishing entries 
Vj = V 0j = {g/m)^N/N B ' (j > 0). 

We denote by C the orthogonal matrix that diagonal- 
izes the problem of coupled oscillators defined in l|A2() . 
such that 



C l C = 1, C\Q 2 + V)C = Q 2 



(A3) 



with a diagonal matrix fl 2 containing the new eigenfrc- 
quencies. Then we can express the old coordinates, mo- 
menta and boson operators in terms of the new normal 
modes: 



/mfi.- . P . 



(A5) 
(A6) 



where d denotes a column vector whose entries are the 
operators dj. In order to evaluate the correlator of the 
cm. mode, we have to relate the old annihilation opera- 
tors dj to the new ones. This is accomplished by inserting 
the relations (|A4|) . 1A5(I into IjAlfl and expressing Q, P 
by d according to l|A6J) . In matrix notation, the result 
reads: 



<mft_ . P , 




(A7) 



Here d' is to be understood as the column- vector contain- 
ing the hermitian conjugate operators as entries. Correla- 
tors of the form (hi(t)b\\ are then obtained by inserting 

the result for do = b± and using the time-evolution of the 

"good" boson operators d(t) = e~ int d, as well as their 
thermal equilibrium expectation values: 




(A8) 



~t 



At T > we would have to consider l^dj(t)djj ^ as 
well. 

We now turn to the actual diagonalization of the prob- 
lem defined in Eq. I)A2|I . i.e. of Q 2 + V with V aj = Vj = 
(g/m) ^N/Nb for j > and all other entries of V equal 
to zero. In fact, for our purposes we only need the spec- 
tral weight of the mode j = of the original problem in 
terms of the new eigenmodes and eigenfrequencies: 



n b 



J=0 



(A9) 



Q = CQ 



(A4) This is related directly to the resolvent 
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R(e) = 



e-(si 2 + v) 



(A10) 



Nb 

R(e - iO) - R(e + iO) = 2 7 ri^G ( 2 /(e - 0|) , (All) 

3=0 

where (ij is the j-th eigenvalue of the matrix SI 2 + V. 
Therefore, we have 



W(u) = — - iO) - R(u 2 + %Q))6(w) , (A12) 
iri 

where we used 5(u> — Clj) — 2lo9(uj)8(ui 2 — ft 2 ). 

The evaluation of the resolvent R(e) is straightforward, 
because the "perturbation" V only connects and j > 0. 
Using G -1 = e — fi 2 , we get 



i?(e) = [G + GUG + GUGUG + . 

oo / N B 

= G 00 V? GjjG 

n=0 \i=l 



J no 



^0 



O I] f _ J fi 2 



(A13) 



3>0 " ~\ 

In the continuum limit, we may use Eq. (0 to evaluate 



V 2 

1 30 



Here the definitions 



/ * ~\ T=0 
e + £<5 — oj 2 



A(e)--r( e )sgn(<5). (A14) 



„ A 




T=0 


2?r— 


(pp) 




m 










dv . 


to J 




e — v 





(A15) 
(A16) 



have been introduced, (where a principal value integral 
is implied in the last line). Using this, we finally obtain 
(with ijAT^ - ljATII) ): 



Appendix B: TWO-PARTICLE GREEN'S 
FUNCTION 



The fermion operators c n in the oscillator are related 
to the auxiliary fermion operators ip(x) via Eq. I|14f) . 
Therefore, the special two-particle Green's function we 
need is given by: 



c„(o)4(i)c^)C(o); 

dxdx'dydy' e *(iV+nV-i*-ny) 



(27T) 



We use the expression of ?/> in terms of the boson field </>, 
Eq. (|lfj|) - to obtain (compare the single-particle expres- 
sions Eqs. HU, ©): 



0)^V, ^(a;, 4)^(^,0) 
1 



(2tt)2 
with the exponent 



e i(w+l)(a;+y-i — » ) £ B(2) 



(B2) 



E, 



(2) 



(/>(y,0)<j>(x',t) 
0(y,O)0(j/',O) 

1 



Here we have introduced 1 



<t>(y,0)<j>(x,t) 

(j>(x' ) t)(t>(x ) t) 

0(x,t)0(y',O) 
(C x + C y + C x , + C y ,) . (B3) 

G, = (^,O) 2 )-(0^,O) 2 ) (o) , 

where the second correlator refers to the non-interacting 
case (without bath) and stems from the factor r (see Eq. 

03). 

Once again, we restrict the further evaluation to the 
special case of T = 0, where some of the expressions 
become slightly simpler. We express the nontrivial cm. 
mode contribution to the correlator of <fi in terms of the 
6^-correlators a(t), 0±(t) introduced above (see Eqs. 
(EH, EH) and (E3J). Since the two-particle Green's func- 
tion involving the oscillator fermion operators c„ is re- 
lated to that of the auxiliary fermion operators ip by 
means of a Fourier integral over the coordinates x, y, x', y' 
(cf. Eq. (|B 1|) ~) . we find it convenient to introduce the 
following abbreviations (similar to our approach for the 
single-particle Green's function): 



w{uj) ~ ^ M p - to 2 - Ap^ + (r(^)/2) 2 



(A17) 



x O = e ix(') yC) = e iy 



(') 



(B4) 



Then we obtain 
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C x = /3+(0) + j9_ (0) - 1 - q(0) (X 2 + X- 2 ), (B5) 

as well as the form of the ^-correlator given already in 
Eq. (EH). 

Therefore, the exponent Em) l|B3|) is given as a sum of 
terms containing a, (3± and binomials constructed out of 
X, X', Y and Y' , as well as non-interacting contributions 

4>(x' , t)<j>(x, 0) ) , that also contain higher powers of 

/(o) 

these variables (see Eq. lUHJO- According to Eq. (|B1|I . 

the two-particle Green's function 

(c n (0)cim(t)ci(0)) (B6) 

is obtained by expanding exp(_E( 2 )) in a power series with 
respect to X, X' , F and Y' , and reading off the coef- 
ficient in front of the component x sl X'- sl 'Y Sn Y'- Sn ' , 
where 51 = 1- (n F + 1) etc. Using JBT}, |E2j and (|B3)l . 
as well as splitting off non-interacting terms, we find that 
the expression to be expanded can be written in the fol- 
lowing form: 

A oe z . (B7) 

Here A yields the correct result for the non-interacting 
case (cf. the application of Wick's theorem, Eq. JHHI), 

r<o v 7 »>o v 7 

rV>0 ^ ' n>0 ^ 7 

and the exponent i? is a rather lengthy expression in- 
volving alt), a(-t) = a*(t), (3±(t), 0±(-t) = 0*±{t) and 
binomials made out of X, X' , Y and Y'; 



In this equation, we have employed the abbreviation 
j3+(t) = P+(t) - e~ luJot . Note that E vanishes if the 
coupling to the bath is switched off. 

Appendix C: LIMITING EXPRESSIONS FOR 
TWO-PARTICLE GREEN'S FUNCTION 

Here we will show how to derive simplified expressions 
for the two-particle Green's function of the previous sec- 
tion, in the limit of weak coupling and for high initial 
excitation. We rewrite Eq. (|B9(I for the weak coupling- 
case, setting a(t) = fi-{t) = 0, /3 + (0) = 1 and keep- 
ing only the (slowly decaying) (3+(t) (compare the main 
text, Sec. I VIII and the corresponding discussion for the 
single-particle case in Sec. II V|) : 

*-*>(f-f) + "«>Gf4)- (C1) 

Here we have set ujq = 0, defining v(t) = [3 + (t)e tuJot = 
e -^t-N-ft/2 _ 1 We wiU algQ get Wo = q inside Aq ^ 

because it can be checked that luq enters the final density 
matrix only in the trivial form exp(iujo(n — n')t). In the 
following, we will shift indices, denoting 51 = I — (np + 1) 
by I, for brevity. Our task then is to obtain the coefficient 
of X l X'- l 'Y n Y'- n ' in the expansion of A expE. 

Let us turn first to the product of the second part of 
Aq (see Eq. (IBSj) ') with the exponential exp(E). As we 
will see, this describes the decay of the excited state, 
while the first part will be connected to heating around 
the Fermi surface. Performing the expansion leads to a 
coefficient 

m 1 +m 2 ,/rai +m 2 
\^ (_-nm 2 +ro 2 ^ ^Q2) 

' mi!r7i2!mi!?7i2! 

where the summation is over all mi 2, ^1,2, n, ft! > 0, 
subject to the constraints: 



E = a(t) (X'Y 



1 



-XY' - 



1 



X'Y' XY 
a*(t)(XY+-±--X'Y- 



a(0) {XX' 



XY 
1 



XX' 



YY> 



X'Y 
1 



X' 



YY' 
Y Y 



/3 + (0)(^ + ^ 



r 



X' 



Mo)( 



Y' X 
X' X 



Y 
Y' 



Y 
X 



Y X' 



-{C x + Cy + C X> + Cy,) . (B9) 



I' 

n 
ri 



n + mi — mi 
n + fhi — mi 
h + rhi + fhi 
n' + mi + mi . 



(C3) 
(C4) 
(C5) 
(C6) 



This requires n — l' = n' — l = m = mi + fhi > 0, i.e. the 
Green's function vanishes unless (I, I') = (n' , n) — m(l, 1) 
for some m > 0, corresponding to the total energy lost to 
the bath (in units of hu>o)- If the given /, n, n' obey this 
relation, then we are left with only three equations for six 
summation indices. The equations < ri = I' + mi — fhi 
and < h' = I + fhi — mi can be fulfilled if fhi < I' + mi 
and mi < I + fhi, respectively. The third equation is 
m 2 + fhi = m. Thus, we have a sum running over m 2 = 
. . . to, toi — . . . I + fhi and fhi — . . . I' + TO2, with 
TO2 = m — mi: 
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(-ir 



mi +m 2 l/ *m 1 +m 2 



(C7) 



This coefficient depends on (m, ?,/'). However, in the 
limit of a high initial excitation, n,n',l,l' 3> 1, we can 
extend the sums over m\ and 777,1 to infinity. This yields: 



Pdocay(ra,t) 



(-1)' 



■{v{t) + v*{t)) m e^ t)+u ' (t) . (C8) 



In this limit, the presently discussed part of the two- 
particle Green's function only depends on the net number 
of quanta m emitted into the bath, and its time-evolution 
is described by the simple expression ilC8(l . For t — > oo, 
we have v(t) — » —1, and thus Pdecay(™,i) - * f^T e_2 - 

The product of the first part of Aowith exp(i?) de- 
scribes how some excitation above the Fermi surface is 
able to influence the states near the surface, by creat- 
ing some "heating" (via the bath-induced interaction be- 
tween particles). Again, expanding leads to a coefficient 
of the same form as Eq. I|C2|I . but with different con- 
straints: 



n 
i 

n 



I + mi — rh 2 
I + fh\ — m% 
n + fhi + fa? 
h + mi + 7712 



(C9) 
(CIO) 
(Cll) 
(C12) 



where I < and h > 0. We find n — V — n' — I 

li / ■ 777 2 + 7772 > 0. The Conditions 1 = 1 + 777 2 — 777! < 

and h = n — [fh\ + 7712) > can be fulfilled if m\ > 
777,2 + I and 777i + 777 2 < 77. Thus, the following sum runs 
over all 7772 = ... 77, 7771 = max(0, V + 1) . . . n — 7772, 

7771 = max(0, 7772 + 1 + 1) . ..I — I' + VH\ + 7772, with 777 2 = 

777i + "72 — m\ +1 — 1': 



, . . | , ,,,,t\mi+?ii2 

fW^ f E M i Vr , (-ir 2+ ^- (CIS) 

' — ' 777l !7772 17771 !777 2 ! 

In the limit 77,77' — * 00, the sums over 7772 and fh\ be- 
come unbounded. Unfortunately, no further simplifica- 
tion is possible. Note that we are interested in small 
I, I', describing states near the Fermi surface, which 
are affected by the high excitation at (77,77'). In this 
limit, the dependence on 77,77' drops out, and we have 

Phoat = PheatG) I 1 1). 



Appendix D: RESOLVENT FOR ARBITRARY 
COUPLING BETWEEN LEVELS 



For the case of damping of all boson modes (transi- 
tions with arbitrary energy transfer, see Sec. IVIIIfl . we 
introduce the resolvent 



WO 



e - (n 2 + V) 



(Dl) 



in terms of which the "spectral weight" W q i q {u>) reads 
(for uj > 0) 



Wg> g (u) = —lmR q , q (cj 2 - i0+) . (D2) 
According to Eq. f77| . the perturbation in the problem 



of coupled oscillators is given by VI 



<i-j 



U, 



3-1 



= vf q q 



{q > 1) j > 1) v = g^N/NB/m), i.e. it always leads to 
transitions between bath modes and boson modes (note 
there are no nonvanishing terms Vjj or V- q .- q >). This al- 
lows us to employ a procedure similar to the previous one 
(Eq. (XT3j|), and to sum the geometric series, in order 
to obtain the following exact result for the resolvent: 



G- q 'G- q q'qf q > f q 



' 2 £,>iG, 



~~ ^2q>l <?$\G-q 

(D3) 

Here Gk = Gkk = (e — fif ) 1 is the unperturbed resol- 
vent. The sum over bath modes j > 1 may be expressed 
in terms of the bath spectrum (see Eq. (J5J in main text). 
This yields, finally: 



Si' 4 fql 



e-(quj ) 2 e- {q'u>o) 2 e- (qu>o) 2 



f 2 a 2 
^ e - {qu ) 2 



q>\ 



(D4) 



We have defined: 



_ 2N 
777 



du> 



I - ~\ T=0 

6 — UJ 2 



(D5) 



The steps up to this point correspond to integrating out 
the bath, which leads to a problem of coupled, damped 
boson modes (whose spectrum is given indirectly, via the 
poles of the resolvent). 

In order to analyze the weak-coupling limit of this ex- 
pression, we treat A as small. Then we obtain the fol- 
lowing result for the off-diagonal terms of the resolvent 
(q' 7^ q), which describe the induced coupling between 
the boson modes: 



Rq' q (e) W Xfq'f q q'q 1— 2 „ 

e- {q'uJo) 2 - \f q ,q' 2 



e - (qcuo) 2 - Xf 2 q 2 



(D6) 
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This expression has been derived by keeping only the 
poles near qtu and q'coo and neglecting terms of order A 2 
in the denominator of Eq. l)D4jl . The prefactor A makes 
these off-diagonal contributions small for small system- 
bath coupling. 

On the other hand, the diagonal contributions may be 
approximated as: 



B„ = — — 



C„ 



D„ 



1 (Y 1 

q \Xw / 



1 fX' 



q \ Y> 

1 fx 

-w 



-it' 



q \Y' 



(E2) 
(E3) 
(E4) 



e-(^ )2-A/ g V 



(D7) 



Here the approximation consisted in dropping poles other 
than that at quo in the sum over q in expression (|D4|1 . 
Inserting this approximation (|D7J| into Eq. (|D2J| for 

W qq , we obtain: 



W qq (uj) sa — Im (uj 2 

1 



i0+ - (qcj ) 2 

— IlTl(5Sl_„ 



\ J-2 2\ 



7r (lu — gt^o — Rei50_ g ) 2 + (Im5f2_ g ) 



Here, A in the first line is to be evaluated ate = 



(D8) 



■i(T 



In the second line, we have furthermore approximated 
W gg (a->) as a Lorentz peak of width Im<5f2_ 9 , where the 
complex frequency shift <5f2_ g of boson mode q is given 
as: 



SSl- g = qfg 



N 



lu(FF) 



muiQ J (<?wo) 2 — «0+ — cj 2 



(D9) 



The Fourier tansform of the weak-coupling result (|D8|) for 
W qq is, therefore, an exponentially decaying oscillation, 
exp[— i(qujQ + Re 5Q_ q )t — Im<5£l_ g |i|]. 



Then we have: 



E {2) ^^A q + B q + C q +D q 
q =i 

^ 1 { [Y\ q [X'^ q 



(E5) 



As explained in Appendix [HI the two-particle Green's 
function 



: n (0)4m(t)cl(0) 



(E6) 



is given by the coefficient of the term 
in the series generated from 



x 0l x 



exp(i?( 2 )), where 51 = 1 — (rip + 1) etc. 

First, the last sum in (|E5(I is split off from £^2) (as it 
does not depend on the time-dependent factors w and u). 
It yields: 



exp 



fc x ,fc„=0 



~x 



(E7) 



Appendix E: EVALUATION OF TWO-PARTICLE 
GREEN'S FUNCTION 



In this appendix, we show the detailed steps needed in 
the actual numerical evaluation of a two-particle Green's 
function. We do this for the case of coupling between 



arbitrary levels (see section IVIII|I , for the Ohmic bath 
at T = 0, evaluated in the weak-coupling approxima- 
tion given by Eq. (|85|l . The resulting approximation 
for the two-particle Green's function can be obtained 
by inserting Eq. (|85|l of the main text into the expres- 
sion for the exponent E( 2 ) (Eq. i|B3|) of Appendix [5J 
where C x = C x i = C y = C y i = in the weak-coupling 
case). In writing down the result, we use the abbrevia- 
tions w = exp(— iu)' t), u = exp(— 7|f|) (and X = exp(ix) 
etc., as before), as well as: 



(El) 



We denote by K(l,l' ,n,h') the coefficient of 
X l X'~ l Y n Y'~ n in the expansion of the remaining 
exponential exp(i?( 2 )) (where that sum has been left out 
in the exponent). According to Eq. 1|E7(I . the desired 
coefficient of the full exponential exp(_E( 2 )) is then given 
by summing the coefficients K over all k x , k y > with 



I = 51 + k x , I = 51 + k x , h = Sn — k y , n — 5n — n, y , 



for given 51, 51' , Sn, 5n' . 

Turning now to the evaluation of K itself, we have 



(E8) 



a: 



q =l \a q =0 q 



(E9) 



and analogous expressions for the other three contribu- 
tions to exp i?( 2 ) , with corresponding sets of exponents 
/3g,7 9 and 5 q . 
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We introduce a = YlqLi 9 a <?' an( ^ analogous definitions 
for /8,7, <5. Then, by looking at Eqs. (|E9|) and the defini- 
tions (|E1|) - (|E4|) . we find that the following relations have 
to be_ fulfilled in order to obtain the desired coefficient K 
of X l X'- V Y h Y'- n ' in the expansion of exp{E (2) ): 

1 = 5-0,1' = a - 7, n = a + /3,h' =-y + 6 (E10) 

Thus, the coefficient K is calculated by adding up 
the following contribution for each set of exponents 
ari, a2, . . . , (3i, . . . , 71, . . . , Si, . . . > that fulfills Eq. 
llETOl) : 

2 l y 

(ll q q a ^) (il^/v) (n 9 9 7 -7 9 !) (n fl «Hi) 

U E, 9 2 (a 9 +/3,+7,+<5q) u; 7+<5-a-/3 > (Ell) 

According to Eq. (|E8(I . we have the upper bounds h — 
a + [3 < Sn and n! = 7 + S < Sn 1 (and, consequently 
from Eq. (|E10(1 . |f|, \l'\ < ra&x(Sn,5n')). Together with 
a q , P q ,~/ q , S q > and the definition of a, this means we 
have to consider only q- values up to q = Sn for a q , (3 q and 
q = Sn' for ~f q , S q ; all exponents for higher q must vanish. 
This justifies the use of the special Ohmic weak-coupling 
form l|E7(l of the exponent, as long as the initial excitation 



(described by Sn, Sn') is not too far from the Fermi level. 
For lower q, we have upper bounds of qa q < Sn (and so 
on). 

In the actual numerical evaluation, the following proce- 
dure is used: Given a maximum excitation of (5n max , a ta- 
ble of coefficients K for all possible values of a, j3, 7, 5 = 
. . . <5rt max is produced, by generating, for each set 
(a,f3,j,S) all possible ai, a 2 , . . . ,/3i, . . . ,71, . . . , Si, . . . 
that fulfill a = Y^Li ^ a q etc - an d adding up the re- 
sulting contributions from l|Ell(l . Since w and u are 
actually functions of time, the contributions are not 
added up immediately but stored in a symbolic form, 
with the numerical prefactor, the w-exponent and the w- 
exponent as entries. The results may finally be used to 
calculate the two-particle Green's function, for any given 
Sl,Sl',5n,5n' < 5n max and any given time t. It is given 
by: 



K(l = Sl+k x , I' = Sl'+k x , h = Sn—ky, h' = Sn'—k y ) , 

kx,k y >0 

(E12) 

where a,(3,j,S are obtained from / etc. by Eq. i|E10|l . 
Note that, because of 1,1' < m&x(Sn, Sn') and h,n' > 0, 
the summation variables k x ,k y are also bounded from 
above. 
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